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Abstract 

(N 

Measurements of flavor ratios of astrophysical neutrino fluxes are sensitive to the two yet un- 
known mixing parameters #13 and S through the combination sin #13 cos 5. We extend previous 
studies by considering the possibility that neutrino fluxes from more than a single type of sources 
will be measured. We point out that, if reactor experiments establish a lower bound on #13, then 
neutrino telescopes might establish an upper bound on | cos 5\ that is smaller than one, and by that 



prove that CP is violated in neutrino oscillations. Such a measurement requires several favorable 
ingredients to occur: (i) #13 is not far below the present upper bound; (ii) The uncertainties in 
012 and 623 are reduced by a factor of about two; (iii) Neutrino fluxes from muon-damped sources 
are identified, and their flavor ratios measured with accuracy of order 10% or better. For the last 
condition to be achieved with the planned km 3 detectors, the neutrino flux should be close to the 
Waxman-Bahcall bound. It motivates neutrino telescopes that are effectively about 10 times larger 
than IceCube for energies of 0(100 TeV), even at the expense of a higher energy threshold. 
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I. INTRODUCTION 

One of the main goals of future neutrino experiments [l| is to observe CP violation in 
neutrino oscillations. The significance of such a measurement goes beyond the determination 
of a fundamental parameter of Nature: it can give further qualitative support to leptogenesis, 
the idea that the observed baryon asymmetry of the Universe has its source in a lepton 
asymmetry generated in neutrino interactions. In some scenarios, it is even quantitatively 
related to leptogenesis. 

Neutrino telescopes [2], such as the IceCube experiment, aim to observe neutrinos coming 
from astrophysical sources. The experiments will provide information on the direction, 
energy, and flavor of the incoming neutrinos. In particular, ratios between fluxes of different 
flavors arriving to the detector can be measured. Ratios between these fluxes at the source 
are predicted by rather robust theoretical considerations. 

The modifications of the flavor ratios between source and detector originate from neutrino 
oscillations. This means that the relations between the fluxes at the source and the fluxes at 
the detector depend on the neutrino parameters in a calculable way. Flavor measurements 
in neutrino telescopes can thus provide information on the neutrino mixing parameters 
. In particular, there is sensitivity to two yet unknown parameters: the 
mixing angle # 13 and the CP violating phase S. 

CP violation in neutrino oscillations can, in principle, be observed via interference terms. 
For neutrinos coming from astrophysical sources, such interference terms are washed out, and 
the measured fluxes are therefore sensitive only to CP conserving parameters. Specifically, 
the measured flavor ratios are sensitive to the combination 

A13 = sin #13 cos 5. (1) 

Since #13 is experimentally bounded from above and known to be small, it is convenient to 
write the flavor ratios in the general form a + &A 13 , where a and b are known functions of 
the two measured parameters, 6*12 and 6*23, but independent of # 13 and 5. The 6A 13 term 
provides a small correction to the zeroth order prediction a. If sin #13 = 0, or if CP violation 
is maximal, i.e. 8 = ir/2 or 3n/2, the correction term is absent. 

If sin #13 is close to the present experimental upper bound, it is likely to be measured in 



near future reactor experiments lOj. In that case, if neutrino telescopes are able to exclude 
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a correction term as large as ±6 sin #13, they will establish that cos 5 7^ ±1 and, by that, will 
discover that CP is violated in neutrino interactions. 

Our goal in this paper is to analyze whether such a discovery of CP violation by neutrino 
telescopes is at all possible. More concretely, we do the following. On the qualitative level, 
we find what types of sources and what types of flavor ratios provide the strongest sensitivity 
to the parameters of interest. On the quantitative level, we estimate the accuracy that is 
required in these measurements and in independent measurements of the mixing angles in 
order to establish that the CP violating phase is different from and from it. Our final 
conclusion is that, with large #13 and near-maximal CP violation, and under some favorable 
circumstances, it may be possible for IceCube (or, more easily, for future, larger detectors) 
to establish CP violation in neutrino interactions. 



II. FLAVOR RATIOS AND MIXING PARAMETERS 



Our goal in this section is to derive analytical expressions for neutrino flavor fluxes that 
can be measured in neutrino telescopes and, in particular, in IceCube. 

Neutrino telescopes can identify the neutrino flavor (a = e, fi, r) via its characteristic 



interaction topology ll|, I12J. IceCube has an energy threshold ~ 100 GeV for detecting 



muon tracks, and ~ 1 TeV for detecting electron- and tau-related showers. Above an energy 
threshold ~ 1 PeV, it is possible to distinguish between the electron-related electromagnetic 
showers and the tau-related hadronic showers. Finally, around E ~ 6.3 PeV, the Glashow 
— ^aaiowtheide^ca^of^eventsaQ. 

We denote the flux of u a + v a measured at the detector by § d a \ the flux of antineutrinos 
v a is denoted by 0^. We consider the following flavor ratios: 

<t>t 

R = # + (2) 

<h d 

s s % < 3 > 

6 d 

T = t (4) 

Below E ~ PeV, only R can be measured. At higher energies, S and perhaps T may become 
available. 
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We denote the flux of u a + v a emitted from the source by <p s a . The relation between (p s a 



and <j)i is given by 



(5) 



where Pp a = P{v a — ^ vp) is the transition probability from a flavor v a at the source to a 
flavor vp at the detector. 

For propagation over astronomical distance scales, the distance-dependent oscillatory 
terms average out, and Pp a depends on mixing parameters only: 



Pf3a — |^ai| 2 |fft| 2 - 



(6) 



Here U is the unitary transformation that relates the neutrino interaction eigenstates v a 
[a = e, [A, r) and mass eigenstates z/j (i = 1, 2, 3): 



Wa) = KM)- 



(7) 



We parametrize the matrix U by three mixing angles, #12, $23 and #13, and three CP violating 
phases, 5, ati and a 2 : 



( 



U 



c 12 c 13 



S12C13 



A 



"S12C23 — Ci2S23Si3e J<5 C12C23 — Si2S22,S\2,e lS S23C13 



V 



S12S23 - ci 2 c 2 3Si 3 e 



-C12S23 - si 2 C2 3 si 3 e C23C13 



ioi/2 



V 



(8) 
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where Cy = cos 9^, s^- = sin^-. It is clear from Eq. §6§ that Pp a is independent of the 
phases «i i2 . It depends on the three mixing angles 6^ and on 5. 

Since it is experimentally known that # 13 is small (see Tabled]), it is convenient to write 
down the flavor transition probabilities to first order in A13 (see Eq. ([!])) Q, 



P ee ~ 1- -sin 2 20i2, 



1 
2 



P eA1 ~ - sin 2 20 i2 cos 2 23 + - sin 20 23 sin 40 i2 A 13 , 



1 

4 



1 ij.fi — 



1 - i (cos 4 023 sin 2 20 12 + sin 2 20 23 ) - - sin 20 23 cos 2 23 sin 40 12 A 13 , 
P er ~ - sin 2 20 12 sin 2 23 - ^ sin 20 23 sin 40 i2 A13, 
PfiT - \ sin 2 2023 (4 - sin 2 20 12 ) + ~ sin40 23 sin40 12 A 13 . 



(9) 



The remaining probabilities can be derived from P a p = Pp a and J2 a P a i3 = 1- 



III. ASTROPHYSICAL NEUTRINO SOURCES AND FLAVOR RATIOS 



We consider two types of sources: 

• "Pion sources" (denoted by sub-index n) provide the following flavor ratios: 

<■>; : o; : <:k 1:2:0. (10) 

As concerns the v e — v e decomposition of <f) e , the situation depends on whether the 
pions are produced mainly by pp or py interactions: 



£ = 1 € 
^2' 0* 




(11) 



"Muon-damped sources" (denoted by sub-index /i) provide the following flavor ratios: 

Of : o; : o: 0:1: 0. (12) 

As concerns the — decomposition of the situation depends on whether the 
pions are produced mainly by pp or p'j interactions: 



1/2 m (13) 
pj. 



The expectation is that all sources where the initial stage of neutrino production is 

rom a "pion" to "muon-damped" flavor de- 



171 ] . If energy losses are mainly due to syn- 



charged pion decays will undergo a transition 
composition at high enough neutrino energies 
chrotron radiation and inverse compton emission, the transition region is expected to span 
about one decade in energy. The actual threshold energy cannot be determined model inde- 
pendently and, furthermore, is likely to differ from source to source. We assume here that, 
nevertheless, the transition is such that it will be possible to separate the neutrino events to 
lower-energy events from pion sources and higher-energy events from muon-damped sources. 

The dependence of the flavor ratios at the detector on the mixing parameters can be 
obtained as follows. One starts from the fluxes at the source (in arbitrary units), Eqs. ( TTOT) . 
(TTTj) . (fl2l) and (TlBl . Then, the fluxes at the detector can be found by using Eq. (]5]) and the 
expressions for the transition probabilities ([6]). Finally, the expressions are put in Eqs. ([2]), 
(H and (H. 



TABLE I: Experimental ranges of mixing angles [1 



Parameter 


Best fit 


lcr range 


'Future' 


sin 2 9 12 


0.31 


0.29-0.33 


0.31 ±0.01 


sin 2 #23 


0.47 


0.40 - 0.55 


0.47 ±0.04 


sin 2 6»i3 


0.00 


< 0.008 


0.022 ± 0.003 



IV. DESCRIPTION OF ANALYSIS 



A. Numerical input 

The current best fit values and la ranges of the mixing angles are given in Tabled] 
By the time that IceCube can carry out the measurements that we discuss in this work, it 
is likely that the knowledge - from other experiments - of the mixing angles will improve. 
Such progress is very significant for our purposes, as we see below. In particular, in order 
that the IceCube measurements will be able, even in principle, to show that 5 ^ 0, it is 
crucial that experiments establish that sin 813 7^ 0. For the sake of our analysis, we assume 
that reactor experiments will measure sin 2 2^13 = 0.090 ± 0.013 [h]]. (This value for 0x3 
corresponds to the current 2a allowed range l].) For 812 and 823 we assume a factor of two 
improvement in the accuracy. The resulting ranges which we use to examine the question 
of whether IceCube can discover CP violation are given in the column labelled 'Future' in 
Table [H As concerns the phase, we assume that it will remain unconstrained. 

To obtain an understanding of the dependence of the flavor ratios on the mixing parameter 
A13, we use the central values for the two measured angles, 812 and 823, and apply the 
approximate relations (Q. We obtain for the pion source 

R n = 0.49 - 0.15Ai3, 
S n = 1.04±0.52A 13 , 

f0.52 + O.28A13 pp. 
(14) 
0.23 + 0.22A 13 pry, 
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and for the muon-damped source 



= 0.62 - 0.49Ai 3 , 
= 0.58 + 0.44A 13 , 

I 0.30 + 0.38A 13 pp, / 
T, = (15) 

[0 pj, 

We emphasize, however, that in our calculations we use the full dependence on the mixing 
angles [see Eq. (JBj)], and not just the leading order (in Ai 3 ) expressions, Eqs. (jSJ), ffl4l) and 

(USD. 



B. Experimental errors 

It is not yet clear whether all of the flavor ratios defined in Section [Til will indeed be 
available at IceCube (or any future neutrino telescope). We assume that R n , R^ and will 
be measured, and consider cases where S n and T M are available or not. 

The goal of this work is not to obtain a detailed realistic estimate of the accuracies that 
are expected in the relevant measurements. Such an estimate depends on both features of 
the astrophysical neutrinos that are not yet known (e.g. the actual total flux), and features 
of the detectors that will only become clear when these neutrinos are observed. The main 
goal here is to find the accuracies that are required in order to establish that CP is violated. 

We thus consider the following experimental accuracies in the measurements of the various 
flavor ratios: 

1. R n : we consider hypothetical accuracies of 5%, 10% or 20%. If the flux is close to the 
Waxman-Bahcall bound, then we expect (9(100) events, and an error of order 10% 
seems realistic; 

2. Sn- In the cases that it is available, we relate the accuracy to that of R^, by 
assuming a Poisson distribution of the number of events for each neutrino flavor. 
We neglect issues of efficiency in detecting tracks versus showers. This leads to 
AS n /S n = y/SJl + S- l ) 2 /(l + R- 1 )(AR 7T /R 7T ). Using central values from Eq. (fll), 
we obtain AS n /S n = 1.2(Ai2 7r /i? 7r ); 
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TABLE II: Scenarios for experimental accuracies 



Scenario 




AR^/R/j, 




(ASJS n ) 


(AT M /T M ) 


(5,5) 


5 


5 


7 


6 


5 


(5,10) 


5 


10 


13 


6 


10 


(5,20) 


5 


20 


27 


6 


20 


(10,10) 


10 


10 


13 


12 


10 


(10,20) 


10 


20 


27 


12 


20 


(20,20) 


20 


20 


27 


24 


20 



3. _R M : we consider hypothetical accuracies which are at best the same as the error on 
R n and at worst 20%; 

4. S^. Following the same line of thought as for S n , we use AS^/S^ = 
^Js^l + S- l ) 2 /(l + R^)(AR IJL /R ll ). Using central values from Eq. (EGO), we obtain 
ASp/S^ = 1.?>(ARJR^ 

5. T M : In the cases that it is available, we assume AT M /T M = AR^/R^. 

The various scenarios can be defined by the assumed accuracies in R v and R^. We 
denote by (a, b) a scenario where the errors are AR^/R^ = a% and AR^/R^ = b%. The six 
scenarios that we consider are presented in Table (TlJ 

We thus consider a hypothetical set of measurements - R, S, T and sin 2 0^ - which 
provide information on 0^ and 5. The statistical procedure by which this information is 
extracted is described in the following section. 



C. Statistical procedure 



E 



Given a measurement of an observable Y r 

,1 2 



(Y) ± oy, we construct x 2 (6'.y,5) 



Y ^ ^ 6> ' J ' < ^ , where 5^(%, 5) represents the theoretical description of the Y observable. 
The uncertainty cry is given in Table [J for sin 2 0y and in Table [TT] for R, S and T. A 



S 



statistical handling of the parameters is performed by analyzing the quantity Ax 2 (%, 5) = 
X 2 - min^ i5 {x 2 }- 

We define the A-dimensional u a% CL acceptance region", for a subset of N out of 
the four mixing parameters (9ij,5), by the region in the N parameter space for which 
^Xmarg < C~ l (a, N). Here Axj^ is obtained by marginalizing A% 2 with respect to the 
4 — N redundant parameters and C~ l {a, N) is the inverse chi-square CDF with N degrees 
of freedom, evaluated at the point a. We have compared this procedure to the more compu- 
tationally demanding FC construction, (as described in 18) and demonstrated, for example, 
in [19]) under the assumption of gaussian measurement errors, for several sample configura- 
tions. We have found a reasonable agreement between our simplified method and the full FC 
routine, with the former tending in general to supply slightly more conservative acceptance 
regions. 

We define the "a% CL acceptance interval", for a specific parameter, by the set of 
parameter values for which the condition Ax^g < C _1 (a, 1) is satisfied, with Ax^g given 
by marginalizing A% 2 with respect to all of the other parameters. 

An "a% CL fraction of coverage" is further defined for a specific parameter as the per- 
centage of the parameter range that is included in the a% CL acceptance interval. The 
lower is this fraction, the stronger is the exclusion power of the experiment with respect to 
the relevant parameter. 

We say that a specific value of a parameter is excluded with a% confidence, if this value is 
not contained in the corresponding a% acceptance interval. This notion will be used below, 
when we discuss the prospects of various measurement scenarios do exclude CP conservation 
in neutrino oscillations. 



V. RESULTS 

A. Neglecting uncertainties in 9\2 and #23 

To understand the abilities and difficulties that are intrinsic to the measurements by 
neutrino telescopes, we first carry out an analysis where 612 and 823 are held fixed at their 
current best fit values. In the next section, we will study the implications of the uncertainties 
in these angles. 
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We begin by choosing specific values for the parameters 9 is and 5, which we call "true 
parameters". Concretely, we assume a true value #13 = 0.15, and consider mainly three 
possibilities for the true value of 5: the two CP conserving ones (S = 0, n) and the maximally 
CP violating one (5 = it/ '2). We evaluate the flux ratios that theoretically correspond to 
these mixing parameters. For the sake of illustration, we assume that the experimental 
measurements will obtain these flux ratios as their central values, with errors as specified for 
each of our six scenarios. We then perform a fit to #13 and 5 (obtaining, of course, the "true 
values" as the best-fit parameters, but with acceptance regions that are different between 
the various scenarios). 

The resulting 90% CL acceptance regions in the # 13 — 5 plane are presented, for the six 
scenarios, in Figs. [TJ [2] and [3J As can be seen in the figures, for some cases, the neutrino 
telescope measurements can mildly improve our knowledge of #13 compared to the reactor 
constraint. 

As concerns S, the 90% CL fraction of coverage in case that all the relevant observables 
will be measured is shown in Fig. HI for true 6*13 = 0.15 and scanning values of true 5 
between and 7r. Since only CP-even quantities are considered, the results for 5 = it + 9 
are equal to those for 5 = 71 — 9. We can make the following statements: 

1. If the neutrino telescope measurements reach the accuracy assumed in this work, they 
are likely to exclude a certain range of 5. 

2. If the Dirac phase is small (that is close to or 7r), the excluded range will be quite 
significant. 

3. The combination of all available observables is usually significantly more efficient than 
partial combinations. 

4. The power of combining measurements is particularly significant as resolutions get 
worse and in the large phase (8 ~ 7r/2) case. 

5. If only i?7r is measured, no range of 5 will be excluded. 

The main question that we are asking is the following: Given a hypothetical situation 
where 5 ~ 7r/2, will IceCube be able to establish CP violation, that is, exclude and ir 
from the acceptance interval in 5? The answer depends of course on which of the various 
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scenarios described in Table HT| if any, will indeed be achieved in the experiment. The main 
lessons that we draw from our calculations are the following: 

(5,5): Measuring and with an accuracy that is significantly better than 10 percent 
will enable a discovery of CP violation in neutrino oscillations. 

(5,10): With this scenario, the sensitivity to CP violation is only marginally affected if 
either or S n are removed from the analysis. Studying the acceptance interval for 5, one 
finds that CP violation may be established even without either T M or S n . This result will be 
further qualified when we elaborate on the scenario, below. 

(10,10): If both T M and S n are measured, with an accuracy ~ 10%, than the required 
accuracy on R v can be somewhat relaxed. 

((5, 10, 20), 20): If the flavor ratios from muon-damped sources cannot be measured with 
an accuracy significantly better than 20%, then even an excellent measurement of flavor 
ratios from pion sources will not exclude CP conservation. 

We learn that the (5, 10) scenario gives a reasonable sense of the minimal required set 
of measurements and accuracies in order that a discovery that CP is violated in neutrino 
oscillations will become possible. Further insight into the role of each of the five observables 
in achieving this goal is given in Fig. [51 depicting the flavor ratios as a function of 5 and 
the x 2 composition for true 5 = n/2. While measurements of and R^ at the assumed 
accuracies suffice to exclude 5 = 7T, at least one of S n or T M needs to be added in order to 
exclude 5 = 0. 

The probability that CP conserving values of 5 will be excluded as a function of the 
true 5, within the four scenarios (5,5), (5,10), (5,20) and (10,10), is shown in Fig. El To 
produce this plot, we generated a large sample (1000) of random sets of observables with 
the prescribed statistics, then checked for each realization whether 5 = or it is contained 
in the resulting acceptance interval. For example, with zero uncertainties in # 12 and 623, the 
conditional probability to exclude CP conservation in the (10,10) scenario given maximal 
phase is about 50%. Note that statistical fluctuations may lead to erroneous exclusion of 
CP conservation even with sin 5 = 0. The fact that the (10,10) scenario is more likely than 
(5,20) to establish CP violation is suggestive for future detector optimizations: If the errors 
on 9 12 and #23 at the time of analysis are significantly reduced, then it may be preferable to 
improve the detection efficiency at the higher range of the spectrum, E > 100 TeV, even at 
the cost of somewhat weaker efficiency at lower energies. 
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B. Taking into account uncertainties in 9 12 and #23 

As a first step in this analysis, we considered the present ranges for 6*12 and 6*23 (see Table 
[I]). The potential of neutrino telescopes to exclude a range of 5 can be seen from Fig. H] 
(upper right panel). The impact of the uncertainties in 9 12 and 9 2 j, can be understood by 
comparing it to the upper left panel. We learn that, with present accuracies, the excluded 
ranges are weaker by 30-50% compared to the idealized case of zero uncertainties. (The 
importance of this ingredient in the analysis was noted in j^.) 

As a second step, we assumed experimental errors on sin 2 # 12 and sin 2 #23 that are reduced 
by a factor of two compared to the present (see Tabled]). The results are shown in Fig. H] 
(lower panel). By comparing to the upper right panel, we learn that such an improvement 
will entail an exclusion power stronger by about 20% compared to the situation that present 
uncertainties remain. 

Concerning the probability that CP violation will be established, we repeat the analysis 
with the present and with the assumed future uncertainties for the four leading scenarios. 
The results are shown in Fig. Without an improvement in the determination of 9\ 2 and #23; 
only the very optimistic scenario (5,5) allows a discovery. With the assumed improvements, 
the more realistic (5,10) scenario also has over 30% probability to make such a discovery. 
The (5,20) and (10,10) scenarios are not powerful enough to do so. 



C. Discussion 



• 0,3, 



A related analysis has been performed previously in Refs. [6|, |8j, which highlighted the 
synergy between neutrino telescopes and terrestrial experiments. The conclusion in Refs. 
6, 8fl regarding the impact of neutrino telescopes on the issue of CP violation is more 

n n 

pessimistic than ours. The main difference lies in the fact that Refs. p, |8( consider the 
information of one type of sources at a time, and indeed we agree with the pessimistic 
conclusion in this case. What we show, however, is that by combining the two types of 
sources that we considered, the ability to exclude CP conservation improves considerably. 
Actually, if this combination of sources is indeed available (and the experimental accuracy is 
similar to or better than our (10,10) scenario), the exclusion power that neutrino telescopes 
have on 5 will be comparable to the proposed superbeams [2l| . (This situation actually 
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reinforces the point made in [6j: since the 5- dependencies of the IceCube and the superbeam 
measurements are different, the information from the two will be complimentary. 
Ref. 



22] points out that variations in the flavor ratios between sources can reach the 
ten percent level and consequently play an important role in the investigation of the mixing 
parameters from astrophysical neutrinos. In particular, the resulting uncertainties may 
wash-out the effects of the Ai 3 terms, especially in the case of low 13 . We agree that 
flavor composition uncertainties at the source would tighten greatly the requirements on 
the experimental precision. There are two reasons, however, why we think that this issue 
may have only limited consequences for our purposes. First, by the time that this analysis 
can be carried out in IceCube, the theoretical analysis of neutrino spectra, which is only 



at its beginning 
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22l | . is likely to improve considerably. In particular, higher quality 
electromagnetic data, from radio to TeV photon energies, will become available. Second, 
our study is relevant only for the case of large #13 where, as we have argued, 10% accuracy 
might be just enough for our purposes if a global analysis of flavor-dependent spectrum will 
be possible. 

The general trends reflected in our results can be simply understood, based on Eqs. ( |T4l) 
and f|T5|) . We rewrite them as follows: 

R n = 0.49 [1 - 0.05(s 13 /0.15) cos 6} , 
S n = 1.04 [1 + 0.08(si 3 /0.15) cos5] , 

= 0.62 [1 - 0.12(s 13 /0.15) cos 5] , 

= 0.58 [1 + 0.11(si 3 /0.15) cos 5] , 
Tjj = 0.30 [1 + 0.19(si 3 /0.15) cos 5] . (16) 

We learn the following: 

• The ratios related to muon-damped sources are more sensitive to the cos (^-dependent 
terms than those related to pion sources; 

• To be sensitive to the cos (^-dependent terms, the accuracy should be of order 10% or 
better; 

• The required accuracy scales with si 3 . If, for example, Si 3 ~ 0.05, sensitivity to cos 5 
will be achieved only with accuracy better than 5%, which seems out of reach for 
IceCube. 
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VI. CONCLUSIONS 



We have studied the potential of combining measurements of flavor ratios in neutrino 
telescopes with observation of 6* 13 7^ by reactor experiments in constraining 5, the CP 
violating phase in the lepton mixing matrix. We reached the following conclusions: 

• Since the neutrino telescopes are sensitive only to the combination A 13 = sin # 13 cos S, 
they can constrain S only if sin # 13 is not too small 6| . 

• Neutrino telescope may exclude at 90% CL up to 30% of the a-priori allowed range 
for 5, even with present accuracies in #12 and 623- 

• Since the Ai 3 -term is maximized in size for cos 5 — ±1, the exclusion region is largest 
if CP is nearly conserved [6j. 

• Reduced uncertainties in #i 2 and # 23 can enlarge the excluded region to about 50% of 
the a-priori allowed range, and give sensitivity even for cos 5 ~ 0. 

• Measuring flavor ratios of fluxes from muon-damped sources will further strengthen 
the exclusion power (compared to measurements based on solely pion sources). Their 
significance is particularly important for cos 5 ~ 0. 

A more specific question that we posed is whether, in case that the CP violating phase S 
is large (~ 7r/2), the measurements of flavor ratios among neutrino fluxes from astrophysical 
sources can establish that the phase is indeed different from or ir, and by that prove that 
CP is violated in neutrino interactions. Our conclusions regarding this question are the 
following: 

• sin 813 must be large, between current 1 — 2a upper bounds. 

• The neutrino flux must not be lower than the Waxman-Bahcall bound. If the flux 
is smaller, a larger neutrino telescope may still achieve this goal, within a reasonable 
time scale 10 years). 

• Neutrino flux from muon-damped sources must be identified, and the related flavor 
ratios measured with accuracy better than 10%. 
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• The uncertainties on 612 and 623 must be reduced by other experiments by a factor of 
about two. 

Even if all these conditions are met, the probability of excluding CP conservation in neutrino 
oscillations is at best 60%. 

The strongest sensitivity to cos 5 arises in flavor ratios related to muon-damped sources. 
On the theoretical side, a more careful study of the transition at high energy from pion-source 
to muon-damped source is important for better understanding of this crucial ingredient in 
our analysis [22]. On Nature's side, the lower the transition energy, and the sharper the tran- 
sition, the higher statistics of events from muon-damped source that will become available 
and, consequently, the better chances are that a neutrino telescope will contribute signifi- 
cantly to understanding CP violation in neutrino oscillations. Finally, on the experimental 
side, a neutrino telescope that is effectively ten times bigger than IceCube, for neutrino 
energy ~ 100 TeV (see Section IV Aft , is well motivated by our arguments. 

The fact that establishing CP violation in IceCube, an experiment under construction, is 
not manifestly impossible is exciting. While a combination of several favorable circumstances 
is required to achieve such a goal, it is worth to refine this analysis, to prepare for a fortunate 
case that these circumstances are fulfilled by the parameters of Nature and by the capabilities 
of neutrino telescopes. 
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FIG. 4: 90% CL fraction of coverage for 5 in the six scenarios denned in Table HH The three panels 
differ in the uncertainties attributed to 6*12 and #23 (see Table H]): (upper left) Zero uncertainties; 
(upper right) Present uncertainties; (bottom) 'Future' uncertainties. 
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It/4 lt/2 3n/4 It 5n/4 3n/2 7jt/4 2jt 

FIG. 5: The (5,10) scenario: (left) The flavor ratios as a function of 5 and their one-sigma range 
(arrows mark the central values corresponding to 5 = 7r/2); (right) The x 2 composition for true 
5 = 7r/2. Both panels correspond to 612 and 623 fixed at their best-fit values, #13 = 0.15. 
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FIG. 6: Probability to exclude CP conservation with 90%CL. The three panels differ in the un- 
certainties attributed to 6\2 and #23 (see Table H]): (upper left) Zero uncertainties; (upper right) 
Present uncertainties; (bottom) 'Future' uncertainties. 
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